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1.  Introduction.  This  paper  discusses  a  new  approach  to  the  feedback  control  of 
aircraft  at  high  angles  of  attack.  This  approach  is  based  on  recent  results  on  control  of 
nonlinear  systems  at  bifurcation  points  [3,  4,  5]  .  The  lateral  dynamics  of  a  slender  air¬ 
craft  studied  by  Ross  [6]  provides  a  convenient  model  for  illustrating  this  application.  It 
is  shown  how  local  bifurcation  control  can  be  used  in  the  stabilization  of  the  so-called 
“wing  rock”  phenomenon.  The  appropriateness  of  this  approach  derives  from  the  fact 
that  unstable  motions  at  high  incidence,  such  as  stall,  spin  entry,  wing  rock  and  roll  cou¬ 
pling,  have  been  attributed  to  bifurcations  in  the  vehicle  dynamic  equations.  The  insta¬ 
bility  can  indeed  be  predicted  via  bifurcation  analysis,  and  depends  critically  on  the  non¬ 
linear  aerodynamic  force  modeling.  Bifurcation  control  results  can  facilitate  the  syn¬ 
thesis  of  nonlinear  stabilizing  feedback  controls  for  aircraft  at  the  onset  of  vortex  separa¬ 
tion;  an  angle  of  attack  limiter  is  not  employed  in  this  approach.  Possible  means  of 
achieving  the  required  control  laws  include  the  control  surface  deflections  and  thrust  vec¬ 
tor  control. 

The  bifurcation  control  technique  considered  here,  developed  in  [3, 4]  ,  provides  an 
analytically-based,  algorithmic  approach  to  the  stabilization  of  vehicles  at  the  onset  of 
flow  separation.  In  classical  terminology,  the  goal  of  the  nonlinear  control  design  is  to 
render  the  stability  boundary  at  high  incidence  “safe”  (as  opposed  to  “dangerous”). 
Practically  the  result  would  be  to  enlarge  the  flight  envelope,  by  ensuring  that  any 
near-stall  sustained  deviation  from  trim  is  stable  and  of  small  amplitude,  and  therefore 
tolerable.  Divergence  cannot  result  from  such  a  state  without  the  introduction  of  a  large 
disturbance.  Thus,  pre-  and  post-stall  dynamics  are  stabilized. 

The  remainder  of  this  paper  is  organized  as  follows.  Section  2  presents  a  brief  dis¬ 
cussion  of  relevant  background  material  on  high  angle  of  attack  dynamic  modeling  and 
associated  nonlinear  effects.  In  Section  3,  the  bifurcation  control  framework  of  [3,  4]  is 
reviewed.  In  Section  4,  application  of  the  bifurcation  control  technique  to  a  model  for 
the  high  incidence  dynamics  of  the  Handley  Page  115  research  aircraft  model  of  Ross  [6] 
is  presented. 

2.  High  Angle  of  Attack  Flight  Dynamics.  Analytical  study  of  rigid  body 
flight  at  high  angle  of  attack  (a)  is  made  difficult  by  essentially  three  main  factors: 
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(i)  Strong  dynamic  coupling  of  longitudinal  and  lateral-directional  motions 
-almost  always  necessitates  the  use  of  a  six  degree- of-freedom  model,  with  state 
variables  representing  angular  orientation  and  rates,  as  well  as  translational 
position  and  velocity; 

(ii)  Nonlinearity  of  the  dynamic  motion  equations.  Since  the  trim  condition 
borders  on  instability  in  the  high  a  regime,  linearized  equations  of  motion  are  of 
little  use  in  predicting  dynamic  behavior;  and 

(iii)  Nonlinear  dependence  of  aerodynamic  forces  on  states.  Vortex  shedding 
and  flow  separation  complicate  aerodynamic  modeling  in  this  flight  regime, 
especially  with  regard  to  wing-body  and  tail-body  interference  effects.  In  the 
supersonic  range,  time-history  effects  also  become  important. 

Although  nonlinear  aerodynamic  modeling  is  difficult  and  highly  configuration- 
dependent,  the  dynamics  of  rigid  body  flight  are  well  understood.  Viewing  an  aircraft  as 
a  rigid  body  of  fixed  mass,  its  motion  is  completely  determined  (once  the  aerodynamic 
forces  are  known)  by  the  following  six  equations,  which  are  expressed  in  the  basic  Euler 
form  with  respect  to  body  axes: 


*  ^  X  ^ 

v  = - -f - 1  vr  -  wq 

m  m 

(2.1) 

•  A  y  ^ 

v  = - -f  —  +  wp  -  %r 

(2.2) 

mm 

•  / *  z 

w  = - +  —  -f  q%  -  vp 

(2-3) 

m  m 

p—  T  (  J  )?r  +  —J — (r  +9<l)  +  —I — (Q  tp)  +  ~f— (ff  r2) 

if  if  if  if  if 

(2.4) 

9 =  r  +  (  r  )pr  -j- (P  -r*)+  ( p  +qr)+-y—(r  pq) 

(2.5) 

r— j  -(  j  )w  +  j  (p  -  qr)  +  j  (p  q*)  + -j-(q  +  tp). 

(2.6) 

Here  u  ,v  ,w  are  the  velocity  components  of  the  center  of  gravity  along  the  x  ,y  and  z 
axes,  respectively;  p  ,q  ,r  are  the  components  of  the  angular  velocity  about  these  axes;  m 
is  the  mass  of  the  aircraft;  /,,/,  and  } ,  are  the  Eulerian  components  of  the  instantane¬ 
ous  gravity  force;  X,Y  and  Z  are  the  total  aerodynamic  forces;  and  L  ,M ,N  are  the 
total  moments  with  respect  to  the  standard  (body)  axes.  The  notation  /, ,  In ,  etc.  refers 
to  the  moments  and  products  of  inertia  relative  to  the  center  of  mass. 

If  the  mass  distribution  of  the  aircraft  is  symmetrical  with  respect  to,  say,  the  xz 
plane,  then  /„=/„=  0  and  the  last  two  terms  in  each  of  the  equations  (2.4)-(2.6)  vanish. 
This  is  a  very  good  approximation  for  most  aircraft  in  steady  flight,  but  may  not  be 
appropriate  if,  for  example,  a  significant  amount  of  fuel  motion  occurs  in  a  tank  (“fuel 
sloshing”). 
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The  aerodynamic  forces  and  moments  acting  on  the  aircraft  are  expressed  in  terms 
of  the  familiar  dimensionless  aerodynamic  coefficients: 

X=  (1/2)  C,pV*S,  Y=  (1/2)  C,pV*S,  Z=  (1/2)  CgpV*Sf  (2.7a) 

L  =  (1/2)  ClPV*Sb  ,  (1/2)  GmpV2Se  ,  N  =  (1/2)  Cu  PV*Sb  ;  (2.7b) 

with  CjyCiy  etc.  being  the  dimensionless  coefficients;  p  the  atmospheric  density;  S  the 
reference  area;  c  the  mean  aerodynamic  chord;  b  the  wing  span;  and  V  the  magnitude 
of  the  velocity  of  the  center  of  gravity.  The  aerodynamic  coefficients  CX,C9)  etc.  above 
are  not  constants  but  depend  on  the  orientation  and  angular  velocity  of  the  vehicle. 
Accurate  analytical  evaluation  of  these  coefficients  is  difficult  and  data  from  wind  tunnel 
tests  is  usually  relied  upon  (when  available).  Note  that  the  effect  of  control  surface 
deflections  is  incorporated  into  the  model  through  the  aerodynamic  coefficients.  The 
control  surface  deflections  correspond  to  the  values  of  the  aileron,  rudder  and  elevator 
angles  S4t  8r  and  6t ,  respectively.  The  thrust  may  also  be  viewed  as  a  control  variable, 
such  as  in  thrust  vector  control. 

Equations  (2.l)-(2.3)  may  be  recast  in  terms  of  the  angle  of  attack  a,  the  sideslip 
angle  /?  and  the  velocity  magnitude  V,  where 

a— tan'1— ,  sin'1  and  V ~(u2  +  v2  +  w2)1^2  .  (2.8) 

This  is  facilitated  by  expressing  the  aerodynamic  forces  X,Y,Z  as 

X  ~L  sin  a  -  D  cospco&a 


Y^Dsinfl 


(2.9) 


Z  =-L  cosct-D  coe/feincr 

where  L  and  D  are  the  net  lift  and  drag  components  of  the  aerodynamic  force: 

L  =  (1/2)  CLP  V*S=  lift 

D—  (1/2)  CDpV^S  =  drag 

The  transformed  equations  (for  a  symmetrical  body)  may  be  found  in  Adams  (7]  . 

It  appears  that  the  first  analytical  study  of  bifurcation  or  jump  phenomena  in  high 
a  flight  dynamics  was  due  to  Schy  and  Hannah  [8]  although  predictions  of  the  utility  of 
the  approach  had  previously  been  made  (Phillips  [9]  ,  Pinsker  [10]  ,  Rhoads  and  Schuler 
[11]  ).  This  analysis  benefited  from  the  report  of  Adams  [7]  which  used  constrained 
minimization  techniques  to  solve  the  nonlinear  algebraic  equations  for  an  equilibrium 
spin.  In  [8]  it  was  shown  how  the  so-called  “pseudosteady-state  solutions”  could  be 
obtained  by  solving  a  single  polynomial  in  roll  rate  with  coefficients  that  were  functions 
of  the  control  inputs.  This  was  followed  by  the  paper  Young,  Schy  and  Johnson  [12]  in 
which  the  inclusion  of  aerodynamics  nonlinear  in  the  angle  of  attack  led  to  the  solution 
of  two  polynomials  in  roll  rate  whose  coefficients  are  functions  of  angle  of  attack  and  the 
control  inputs.  Analytical  results  were  compared  with  calculated  time  histories  to 
demonstrate  the  validity  of  the  method  for  predicting  jump-like  instabilities.  This  work 


6 


Only  the  linear,  quadratic  and  cubic  terms  in  an  applied  feedback  «(*)  have  poten¬ 
tial  for  influencing  the  stability.  Stability  may  be  checked  by  computing  so-called  bifur¬ 
cation  coefficients.  Public  domain  software  packages  exist  which  perform  the  bifurca~ 
tion  analysis,  as  well  as  follow  the  bifurcation  branches  in  one  or  two  parameters.  We 
mention,  in  particular,  the  BACTM  software  of  Carroll  and  Mehra  [13]  ,  the  package 
BIFOR2  of  B.  Hassard  [16]  ,  SUNY  Buffalo,  and  the  interactive  package  AUTO  [17] 
authored  by  E.  Doedel  at  Concordia  University,  Montreal.  Bifurcation  analysis  software 
can  be  used,  in  conjunction  with  simulation  tools,  to  obtain  a  global  picture  of  the 
efficacy  of  a  particular  control  law  design. 

Theorem.  Let  the  uncontrolled  system  ((3.3)  with  v  =  0)  undergo  an  unstable  Hopf 
bifurcation  to  periodic  solutions .  Then  there  is  a  smooth  feedback  u(x)  with  u  (0)  =  0 
which  solves  the  local  Hopf  bifurcation  stabilization  problem  for  Eq.  (3.3)  provided  that 

o  f  Re  {  -  2lQ0(r  ,  jLfr)  +  IQ0(F,  i  (2  iut  I  -  L0i )"‘7) 

+  jl[2Lir  +  hr]  }.  (3.7) 

A  whole  family  of  conditions  similar  to  Eq.  (3.7)  may  be  derived,  each  corresponding 
to  a  particular  family  of  nonlinear  stabilizing  control  laws  u  (x ).  The  results  are  available 
for  both  static  and  Hopf  bifurcations  [18]  .  Eq.  (3.7)  indicates  the  relative  computational 
ease  with  which  the  results  may  be  applied;  only  a  Taylor  series  expansion  of  the 
dynamic  equations  and  standard  spectral  calculations  are  needed. 

4.  Wing  Rock  Stabilization  for  the  HP115.  Nonlinear  feedback  control  laws 
ensuring  a  “safe”  stability  boundary  at  high  angle  of  attack  may  be  designed  using  the 
bifurcation  control  results  discussed  in  Section  3.  Use  of  this  approach  has  several 
advantages  over  most  standard  approaches  to  high  angle  of  attack  flight  control:  (i) 
Rather  than  avoiding  stall  through  angle  of  attack  limiters,  the  aim  is  to  use  feedback  to 
permit  stable  flight  in 'the  stall  and  post-stall  regime;  (ii)  All  nonlinear  stabilizing  feed¬ 
back  controls  can  be  parametrized,  allowing  consideration  of  other  criteria  besides  local 
asymptotic  stability,  such  as  size  of  the  domain  of  stability.  (The  theory  provides  a  way 
to  estimate  the  achieved  domain  of  stability.);  and  (iii)  Stall  stabilization  may  well  be 
preferable  to  deferring  stall  to  higher  angles  of  attack,  since,  as  noted  by  Babister  [2]  , 
the  latter  often  leads  to  a  much  more  violent  stall  at  a  higher  angle  of  attack. 

To  indicate  the  type  of  aircraft  control  problem  which  can  benefit  directly  from  the 
bifurcation  control  approach,  we  mention  thrust  vector  control  as  a  means  for  stabilizing 
aircraft  longitudinal  and/or  lateral  dynamics  at  high  incidence.  We  note  that  in  a  recent 
paper  [19]  ,  Ashley  has  concluded  that  thrust  vector  control  will  be  crucial  in  subsonic 
maneuvering  at  extreme  angles  of  attack. 

Consider  the  following  model  [6,20]  for  the  lateral  dynamics  of  the  HP115  aircraft 
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laws  for  a  control  system 

x=f„(x,u).  (3.3) 

In  (3.3),  v  is  a  control  function  and  p  is,  as  before,  a  bifurcation  parameter.  The  main 
assumption  of  the  approach  is  simply  that  the  ‘uncontrolled’  system  (obtained  by  setting 
u  =  0  in  (3.3)) 

*  — /  m(x  >0)  (3-4) 

undergoes  a  bifurcation  as  p  is  varied.  The  most  important  benefit  of  the  approach  is  an 
algorithmic  design  aid  which  generates,  in  a  computationally  efficient  way,  all  the  purely 
nonlinear  feedback  control  laws  u  =  u(r)  ensuring  stability  of  the  new  (“bifurcated”) 
solutions  for  the  controlled  system.  In  the  case  of  a  Hopf  bifurcation,  for  instance,  one 
can  design  feedbacks  ensuring  that  for  any  value  of  the  parameter  p  in  some  domain 
around  the  critical  value,  the  state  converges  either  to  the  original  equilibrium  solution 
or  to  a  nearby  (stable)  periodic  solution  of  small  (nearly  -  zero)  amplitude.  Since  the 
amplitude  of  the  bifurcated  periodic  orbits  is  small  for  small  p,  this  type  of  (stable)  oscil¬ 
latory  behavior  may  be  tolerable  from  a  practical  point  of  view,  and  is  certainly  prefer¬ 
able  to  divergence  or  uncontrolled  wing  rock. 

Our  approach  to  the  local  feedback  stabilization  problem  has  the  novel  feature  that 
it  facilitates  the  derivation  of  generally  valid  analytical  criteria  for  stabilizability,  as  well 
as  specific  stabilizing  feedback  controls.  This  is  possible  through  use  of  bifurcation  for¬ 
mulae  which  involve  only  Taylor  series  expansion  of  the  vector  field  and  eigenvector 
computations.  To  illustrate  the  generality  of  the  approach  and  the  simplicity  of  the  asso¬ 
ciated  calculations,  a  representative  result  is  given  next. 

Consider  a  one- parameter  family  of  nonlinear  control  systems  (3.3)  where  *  €  R* ,  * 
is  a  scalar  control,  p  is  a  real-valued  parameter,  and  the  vector  field  /  „  is  sufficiently 
smooth.  Suppose  that  for  tsO,  Eq.  (3.3)  has  an  equilibrium  point  x0 (p)  which  depends 
smoothly  on  p.  The  linearization  of  Eq.  (3.3)  at*  =  0,  «  =  p  =  0  is  given  by 

x  —  Ax  +  bu  (3.5) 

where  A  -^-(0,0)  and  b  :=  -^-(0,0).  The  nature  of  the  stabilizing  nonlinear  feedack 
dx  du 

controls  given  by  the  theory  depends  on  the  controllability  properties  of  this  linear  sys¬ 
tem.  Rewrite  Eq.  (3.3)  in  the  series  form 

x  =  L0x  +  «7  +  vLji  +  Qo{x,x) 

+  <70(z  ,x  ,x )  +  *  *  •  (3-6) 

where  the  terms  not  written  explicitly  are  of  higher  order  in  x,  v  and  p  than  those 
which  are.  Thus  L0  and  L,  are  square  matrices,  7  is  a  constant  vector,  Q d(x,z)  is  a  qua¬ 
dratic  form  generated  by  a  symmetric  bilinear  form  Q o{x,v)  giving  the  second  order  (in 
x )  terms  at  u  =  0,  p  =  0,  and  C0[x  ,z  ,* )  is  a  cubic  form  generated  by  a  symmetric  tri- 
linear  form  CJ(x,y  ,z)  giving  the  third  order  (in  x)  terms  at  «  =  0,  p  =  0.  (Note  that  L0 
is  simply  A  of  Eq.  (3.5),  and  7  corresponds  to  b .)  Denote  by  r  the  right  (column)  and  by 
I  the  left  (row)  eigenvector  of  L0  with  eigenvalue  iue.  Normalize  by  setting  the  first 
component  of  r  to  1  and  then  choose  /  so  that  lr  =  1. 


10 


[17]  E.  Docdel,  AUTO  User’s  Manual,  Computer  Science  Dept.,  Concordia  Univ., 

Montreal,  May  1984. 

[18]  J.-H.  Fu,  Bifurcation  Control  and  Feedback  Stabilization,  Ph.D.  Dissertation, 

Dept,  of  Elec.  Eng.,  Univ.  of  Maryland,  College  Park,  1988. 

[19]  H.  Ashley,  “On  the  feasibility  of  low-speed  fighter  maneuvers  involving  extreme 

angles  of  attack,”  Tech.  Rept.,  Dept,  of  Aeronautics  and  Astronautics,  Stan¬ 
ford  Univ.,  to  be  published. 

[20]  J.E.  Cochran,  Jr.  and  C.-S.  Ho,  “Stability  of  aircraft  motion  in  critical  cases,”  J. 

Guidance  0,  pp.  272-279,  1983. 


7 


at  high  incidence: 

*a  a  COSO 

p  =  srna  p  -  cosa  r  +  — 
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+  -£<*10  +  «/  +  n,JL-  +  nr-pr)} 


.  PSV2sI,  ' 

r  ~  "j~j  _  j  2  {  n  +  n  +  nP  "y  +  nr 


ra 

V 


+  ~U,f  +  is  + 


(4.1a) 


(4.1b) 


(4.1c) 


9  =  P  +  tana  r  (4. Id) 

Here,  a  is  the  (assumed  constant)  angle  of  attack,  0  and  4>  are  the  angles  of  sideslip  and 
roll,  respectively,  and  p  and  r  are  the  rates  of  roll  and  yaw,  respectively.  The  reader  is 
referred  to  [6,  20]  for  definitions  of  the  parameters  appearing  in  (4.1),  though  we  note 
that  aerodynamic  coefficients  cubic  in  the  sideslip  angle  0  appear,  with  various  parame¬ 
ters  depending  as  well  on  the  angle  of  attack  a.  Precise  forms  of  these  dependencies  are 
given  by  Ross  [6]  ,  which  are  derived  by  approximating  flight  test  data.  Specifically,  for 
a  particular  configuration  and  flight  condition,  we  have 


yr  =  0.014  +  0.505a  -  0.47a2 

(4.2a) 

/,  =-0.132  +  0.08a 

(4.2b) 

np  =  0.0125a  -  0.938a2 

(4.2c) 

II 

O 

(4.2d) 

lr  =0.006  +  0.54a 

(4.2e) 

nr  =  -0.351  -  0.089a 

(4.2f) 

Other  (nonaerodynamic)  parameters  are  set  as  follows  [6] 
V  =  71.25  m/sec,  m  =  2154  kg,  a  —  3.05  m,  S  =  40.18  m2, 

I«  =  25430  kg  m2,  I„  =  1615  kg  m2. 

:  p  =  0.906  kg/ms, 
I„  =  2182  kg  m2, 

The  parameters  which  remain  to  be  specified  in  Eq.  (4.1)  are  the  linear  and  cubic 
coefficients  in  the  aerodynamic  moments  in  sideslip,  roll  and  yaw.  These  are  yj,  y8, 

/s,  and  n,,  ns,  respectively.  The  Hopf  bifurcation  analysis  package  BIFOR2  [16]  was 
used  to  numerically  investigate  the  influence  of  these  parameters  on  the  stability  of  a 
Hopf  bifurcation  in  the  dynamics.  The  results  showed  a  strong  dependence  of  the 
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stability  coefficient  fa  on  the  nonlinear  yaw  parameter  ns.  The  values  of  these  parame¬ 
ters  obtained  from  [6,  20]  are:  yj  =  -0.191,  ys  =  -1.958,  lt  =  -0.184,  /3  =  -0.0,  n,  =  0.07, 
ns  =  2.85.  In  particular,  we  note  that  from  Figure  6  of  [6]  ,  the  values  of  n,  and  na  are 
positive  for  the  nominal  configuration.  Thus,  the  negative  values  for  these  parameters 
quoted  in  [20]  do  not  appear  to  agree  with  the  actual  data.  It  is  therefore  not  surprising 
that  the  calculations  we  have  performed  do  not  coincide  exactly  to  those  of  [20]  .  Note 
also  that  we  have  followed  [6]  in  assuming  an  air  density  value  of  p  =  0.906,  while  a 
value  of  0.5  is  quoted  in  [20]  . 

For  the  nominal  parameter  values  given  in  the  foregoing,  the  program  BIFOR2  com¬ 
putes  a  Hopf  bifurcation  from  the  origin  for  Eq.  (4.1)  with  a  critical  angle  of  attack 
ac  =  0.295  rad  =  16.9  deg,  in  basic  agreement  with  the  references.  The  bifurcation  is 
supercritical,  with  the  bifurcation  stability- coefficient  fa  —  -4.1-  The  calculations  for  the 
nominal  system  also  reveal  a  difficulty  with  the  assumed  aerodynamic  model:  the  eigen¬ 
values  cross  the  imaginary  axis  from  the  right  half  of  the  complex  plane  into  the  left  half 
plane  as  or  is  increased  through  its  critical  value.  For  the  nominal  values  above,  BIFOR2 
determines  that  the  rate  of  change  of  the  real  part  of  these  eigenvalues  with  respect  to  a 
is  -0.937.  The  noncritical  eigenvalues  of  (4.1)  at  the  bifurcation  point  are  found  to  be 
-0.64  and  -0.05,  which  compare  well  with  [20]  when  the  differences  in  the  parameter 
values  used  are  taken  into  consideration. 

In  order  to  illustrate  application  of  the  bifurcation  control  results  discussed  in  Sec¬ 
tion  3,  we  consider  the  effect  of  the  parameter  ns  on  stability  of  the  Hopf  bifurcation. 
Using  a  formula  for  fa  given  in  [3] ,  we  have  determined  that  fa  for  this  problem  depends 
linearly  on  the  coefficient  ns.  This  fact  is  easily  obtained  since  only  cubic  non  linearities 
in  one  component  of  the  state  (0)  appear  in  (4.1).  In  particular,  we  find  that  increasing 
the  value  of  n3  from  2.65  results  in  a  more  negative  stability  coefficient  fa,  and  hence  in 
better  local  stability  properties.  In  contrast,  decreasing  nt  increases  fa,  possibly  leading 
to  a  subcritical  (unstable)  bifurcation.  For  example,  taking  ns  =  -10  results  in  fa  =  2.56, 
and  hence  an  unstable  bifurcation. 
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